3.21 \(\int \frac {\sin ^6(x)}{(a+a \sin (x))^3} \, dx\)

Optimal. Leaf size=101 \[ -\frac {23 x}{2 a^3}+\frac {136 \cos ^3(x)}{15 a^3}-\frac {136 \cos (x)}{5 a^3}+\frac {23 \sin ^3(x) \cos (x)}{3 \left (a^3 \sin (x)+a^3\right )}+\frac {23 \sin (x) \cos (x)}{2 a^3}+\frac {\sin ^5(x) \cos (x)}{5 (a \sin (x)+a)^3}+\frac {13 \sin ^4(x) \cos (x)}{15 a (a \sin (x)+a)^2} \]

[Out]

-23/2*x/a^3-136/5*cos(x)/a^3+136/15*cos(x)^3/a^3+23/2*cos(x)*sin(x)/a^3+1/5*cos(x)*sin(x)^5/(a+a*sin(x))^3+13/
15*cos(x)*sin(x)^4/a/(a+a*sin(x))^2+23/3*cos(x)*sin(x)^3/(a^3+a^3*sin(x))

________________________________________________________________________________________

Rubi [A]  time = 0.23, antiderivative size = 101, normalized size of antiderivative = 1.00, number of steps used = 8, number of rules used = 6, integrand size = 13, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.462, Rules used = {2765, 2977, 2748, 2635, 8, 2633} \[ -\frac {23 x}{2 a^3}+\frac {136 \cos ^3(x)}{15 a^3}-\frac {136 \cos (x)}{5 a^3}+\frac {23 \sin ^3(x) \cos (x)}{3 \left (a^3 \sin (x)+a^3\right )}+\frac {23 \sin (x) \cos (x)}{2 a^3}+\frac {\sin ^5(x) \cos (x)}{5 (a \sin (x)+a)^3}+\frac {13 \sin ^4(x) \cos (x)}{15 a (a \sin (x)+a)^2} \]

Antiderivative was successfully verified.

[In]

Int[Sin[x]^6/(a + a*Sin[x])^3,x]

[Out]

(-23*x)/(2*a^3) - (136*Cos[x])/(5*a^3) + (136*Cos[x]^3)/(15*a^3) + (23*Cos[x]*Sin[x])/(2*a^3) + (Cos[x]*Sin[x]
^5)/(5*(a + a*Sin[x])^3) + (13*Cos[x]*Sin[x]^4)/(15*a*(a + a*Sin[x])^2) + (23*Cos[x]*Sin[x]^3)/(3*(a^3 + a^3*S
in[x]))

Rule 8

Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]

Rule 2633

Int[sin[(c_.) + (d_.)*(x_)]^(n_), x_Symbol] :> -Dist[d^(-1), Subst[Int[Expand[(1 - x^2)^((n - 1)/2), x], x], x
, Cos[c + d*x]], x] /; FreeQ[{c, d}, x] && IGtQ[(n - 1)/2, 0]

Rule 2635

Int[((b_.)*sin[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> -Simp[(b*Cos[c + d*x]*(b*Sin[c + d*x])^(n - 1))/(d*n),
x] + Dist[(b^2*(n - 1))/n, Int[(b*Sin[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n, 1] && Integer
Q[2*n]

Rule 2748

Int[((b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_) + (d_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> Dist[c, Int[(b*S
in[e + f*x])^m, x], x] + Dist[d/b, Int[(b*Sin[e + f*x])^(m + 1), x], x] /; FreeQ[{b, c, d, e, f, m}, x]

Rule 2765

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Sim
p[((b*c - a*d)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*(c + d*Sin[e + f*x])^(n - 1))/(a*f*(2*m + 1)), x] + Dist[1/
(a*b*(2*m + 1)), Int[(a + b*Sin[e + f*x])^(m + 1)*(c + d*Sin[e + f*x])^(n - 2)*Simp[b*(c^2*(m + 1) + d^2*(n -
1)) + a*c*d*(m - n + 1) + d*(a*d*(m - n + 1) + b*c*(m + n))*Sin[e + f*x], x], x], x] /; FreeQ[{a, b, c, d, e,
f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0] && LtQ[m, -1] && GtQ[n, 1] && (IntegersQ
[2*m, 2*n] || (IntegerQ[m] && EqQ[c, 0]))

Rule 2977

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])*((c_.) + (d_.)*sin[(e_
.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp[((A*b - a*B)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*(c + d*Sin[e + f*x]
)^n)/(a*f*(2*m + 1)), x] - Dist[1/(a*b*(2*m + 1)), Int[(a + b*Sin[e + f*x])^(m + 1)*(c + d*Sin[e + f*x])^(n -
1)*Simp[A*(a*d*n - b*c*(m + 1)) - B*(a*c*m + b*d*n) - d*(a*B*(m - n) + A*b*(m + n + 1))*Sin[e + f*x], x], x],
x] /; FreeQ[{a, b, c, d, e, f, A, B}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0] && LtQ
[m, -2^(-1)] && GtQ[n, 0] && IntegerQ[2*m] && (IntegerQ[2*n] || EqQ[c, 0])

Rubi steps

\begin {align*} \int \frac {\sin ^6(x)}{(a+a \sin (x))^3} \, dx &=\frac {\cos (x) \sin ^5(x)}{5 (a+a \sin (x))^3}-\frac {\int \frac {\sin ^4(x) (5 a-8 a \sin (x))}{(a+a \sin (x))^2} \, dx}{5 a^2}\\ &=\frac {\cos (x) \sin ^5(x)}{5 (a+a \sin (x))^3}+\frac {13 \cos (x) \sin ^4(x)}{15 a (a+a \sin (x))^2}-\frac {\int \frac {\sin ^3(x) \left (52 a^2-63 a^2 \sin (x)\right )}{a+a \sin (x)} \, dx}{15 a^4}\\ &=\frac {\cos (x) \sin ^5(x)}{5 (a+a \sin (x))^3}+\frac {13 \cos (x) \sin ^4(x)}{15 a (a+a \sin (x))^2}+\frac {23 \cos (x) \sin ^3(x)}{3 \left (a^3+a^3 \sin (x)\right )}-\frac {\int \sin ^2(x) \left (345 a^3-408 a^3 \sin (x)\right ) \, dx}{15 a^6}\\ &=\frac {\cos (x) \sin ^5(x)}{5 (a+a \sin (x))^3}+\frac {13 \cos (x) \sin ^4(x)}{15 a (a+a \sin (x))^2}+\frac {23 \cos (x) \sin ^3(x)}{3 \left (a^3+a^3 \sin (x)\right )}-\frac {23 \int \sin ^2(x) \, dx}{a^3}+\frac {136 \int \sin ^3(x) \, dx}{5 a^3}\\ &=\frac {23 \cos (x) \sin (x)}{2 a^3}+\frac {\cos (x) \sin ^5(x)}{5 (a+a \sin (x))^3}+\frac {13 \cos (x) \sin ^4(x)}{15 a (a+a \sin (x))^2}+\frac {23 \cos (x) \sin ^3(x)}{3 \left (a^3+a^3 \sin (x)\right )}-\frac {23 \int 1 \, dx}{2 a^3}-\frac {136 \operatorname {Subst}\left (\int \left (1-x^2\right ) \, dx,x,\cos (x)\right )}{5 a^3}\\ &=-\frac {23 x}{2 a^3}-\frac {136 \cos (x)}{5 a^3}+\frac {136 \cos ^3(x)}{15 a^3}+\frac {23 \cos (x) \sin (x)}{2 a^3}+\frac {\cos (x) \sin ^5(x)}{5 (a+a \sin (x))^3}+\frac {13 \cos (x) \sin ^4(x)}{15 a (a+a \sin (x))^2}+\frac {23 \cos (x) \sin ^3(x)}{3 \left (a^3+a^3 \sin (x)\right )}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]  time = 0.10, size = 191, normalized size = 1.89 \[ \frac {\left (\sin \left (\frac {x}{2}\right )+\cos \left (\frac {x}{2}\right )\right ) \left (24 \sin \left (\frac {x}{2}\right )-690 x \left (\sin \left (\frac {x}{2}\right )+\cos \left (\frac {x}{2}\right )\right )^5-405 \cos (x) \left (\sin \left (\frac {x}{2}\right )+\cos \left (\frac {x}{2}\right )\right )^5+5 \cos (3 x) \left (\sin \left (\frac {x}{2}\right )+\cos \left (\frac {x}{2}\right )\right )^5+45 \sin (2 x) \left (\sin \left (\frac {x}{2}\right )+\cos \left (\frac {x}{2}\right )\right )^5+1576 \sin \left (\frac {x}{2}\right ) \left (\sin \left (\frac {x}{2}\right )+\cos \left (\frac {x}{2}\right )\right )^4+112 \left (\sin \left (\frac {x}{2}\right )+\cos \left (\frac {x}{2}\right )\right )^3-224 \sin \left (\frac {x}{2}\right ) \left (\sin \left (\frac {x}{2}\right )+\cos \left (\frac {x}{2}\right )\right )^2-12 \left (\sin \left (\frac {x}{2}\right )+\cos \left (\frac {x}{2}\right )\right )\right )}{60 (a \sin (x)+a)^3} \]

Antiderivative was successfully verified.

[In]

Integrate[Sin[x]^6/(a + a*Sin[x])^3,x]

[Out]

((Cos[x/2] + Sin[x/2])*(24*Sin[x/2] - 12*(Cos[x/2] + Sin[x/2]) - 224*Sin[x/2]*(Cos[x/2] + Sin[x/2])^2 + 112*(C
os[x/2] + Sin[x/2])^3 + 1576*Sin[x/2]*(Cos[x/2] + Sin[x/2])^4 - 690*x*(Cos[x/2] + Sin[x/2])^5 - 405*Cos[x]*(Co
s[x/2] + Sin[x/2])^5 + 5*Cos[3*x]*(Cos[x/2] + Sin[x/2])^5 + 45*(Cos[x/2] + Sin[x/2])^5*Sin[2*x]))/(60*(a + a*S
in[x])^3)

________________________________________________________________________________________

fricas [A]  time = 0.51, size = 158, normalized size = 1.56 \[ \frac {10 \, \cos \relax (x)^{6} - 15 \, \cos \relax (x)^{5} - {\left (345 \, x + 839\right )} \cos \relax (x)^{3} - 140 \, \cos \relax (x)^{4} - {\left (1035 \, x - 668\right )} \cos \relax (x)^{2} + 6 \, {\left (115 \, x + 233\right )} \cos \relax (x) + {\left (10 \, \cos \relax (x)^{5} + 25 \, \cos \relax (x)^{4} - {\left (345 \, x - 724\right )} \cos \relax (x)^{2} - 115 \, \cos \relax (x)^{3} + 6 \, {\left (115 \, x + 232\right )} \cos \relax (x) + 1380 \, x - 6\right )} \sin \relax (x) + 1380 \, x + 6}{30 \, {\left (a^{3} \cos \relax (x)^{3} + 3 \, a^{3} \cos \relax (x)^{2} - 2 \, a^{3} \cos \relax (x) - 4 \, a^{3} + {\left (a^{3} \cos \relax (x)^{2} - 2 \, a^{3} \cos \relax (x) - 4 \, a^{3}\right )} \sin \relax (x)\right )}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(x)^6/(a+a*sin(x))^3,x, algorithm="fricas")

[Out]

1/30*(10*cos(x)^6 - 15*cos(x)^5 - (345*x + 839)*cos(x)^3 - 140*cos(x)^4 - (1035*x - 668)*cos(x)^2 + 6*(115*x +
 233)*cos(x) + (10*cos(x)^5 + 25*cos(x)^4 - (345*x - 724)*cos(x)^2 - 115*cos(x)^3 + 6*(115*x + 232)*cos(x) + 1
380*x - 6)*sin(x) + 1380*x + 6)/(a^3*cos(x)^3 + 3*a^3*cos(x)^2 - 2*a^3*cos(x) - 4*a^3 + (a^3*cos(x)^2 - 2*a^3*
cos(x) - 4*a^3)*sin(x))

________________________________________________________________________________________

giac [A]  time = 0.97, size = 99, normalized size = 0.98 \[ -\frac {23 \, x}{2 \, a^{3}} - \frac {9 \, \tan \left (\frac {1}{2} \, x\right )^{5} + 36 \, \tan \left (\frac {1}{2} \, x\right )^{4} + 84 \, \tan \left (\frac {1}{2} \, x\right )^{2} - 9 \, \tan \left (\frac {1}{2} \, x\right ) + 40}{3 \, {\left (\tan \left (\frac {1}{2} \, x\right )^{2} + 1\right )}^{3} a^{3}} - \frac {4 \, {\left (75 \, \tan \left (\frac {1}{2} \, x\right )^{4} + 330 \, \tan \left (\frac {1}{2} \, x\right )^{3} + 530 \, \tan \left (\frac {1}{2} \, x\right )^{2} + 355 \, \tan \left (\frac {1}{2} \, x\right ) + 86\right )}}{15 \, a^{3} {\left (\tan \left (\frac {1}{2} \, x\right ) + 1\right )}^{5}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(x)^6/(a+a*sin(x))^3,x, algorithm="giac")

[Out]

-23/2*x/a^3 - 1/3*(9*tan(1/2*x)^5 + 36*tan(1/2*x)^4 + 84*tan(1/2*x)^2 - 9*tan(1/2*x) + 40)/((tan(1/2*x)^2 + 1)
^3*a^3) - 4/15*(75*tan(1/2*x)^4 + 330*tan(1/2*x)^3 + 530*tan(1/2*x)^2 + 355*tan(1/2*x) + 86)/(a^3*(tan(1/2*x)
+ 1)^5)

________________________________________________________________________________________

maple [A]  time = 0.10, size = 174, normalized size = 1.72 \[ -\frac {3 \left (\tan ^{5}\left (\frac {x}{2}\right )\right )}{a^{3} \left (\tan ^{2}\left (\frac {x}{2}\right )+1\right )^{3}}-\frac {12 \left (\tan ^{4}\left (\frac {x}{2}\right )\right )}{a^{3} \left (\tan ^{2}\left (\frac {x}{2}\right )+1\right )^{3}}-\frac {28 \left (\tan ^{2}\left (\frac {x}{2}\right )\right )}{a^{3} \left (\tan ^{2}\left (\frac {x}{2}\right )+1\right )^{3}}+\frac {3 \tan \left (\frac {x}{2}\right )}{a^{3} \left (\tan ^{2}\left (\frac {x}{2}\right )+1\right )^{3}}-\frac {40}{3 a^{3} \left (\tan ^{2}\left (\frac {x}{2}\right )+1\right )^{3}}-\frac {23 \arctan \left (\tan \left (\frac {x}{2}\right )\right )}{a^{3}}-\frac {8}{5 a^{3} \left (\tan \left (\frac {x}{2}\right )+1\right )^{5}}+\frac {4}{a^{3} \left (\tan \left (\frac {x}{2}\right )+1\right )^{4}}+\frac {8}{3 a^{3} \left (\tan \left (\frac {x}{2}\right )+1\right )^{3}}-\frac {8}{a^{3} \left (\tan \left (\frac {x}{2}\right )+1\right )^{2}}-\frac {20}{a^{3} \left (\tan \left (\frac {x}{2}\right )+1\right )} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(sin(x)^6/(a+a*sin(x))^3,x)

[Out]

-3/a^3/(tan(1/2*x)^2+1)^3*tan(1/2*x)^5-12/a^3/(tan(1/2*x)^2+1)^3*tan(1/2*x)^4-28/a^3/(tan(1/2*x)^2+1)^3*tan(1/
2*x)^2+3/a^3/(tan(1/2*x)^2+1)^3*tan(1/2*x)-40/3/a^3/(tan(1/2*x)^2+1)^3-23/a^3*arctan(tan(1/2*x))-8/5/a^3/(tan(
1/2*x)+1)^5+4/a^3/(tan(1/2*x)+1)^4+8/3/a^3/(tan(1/2*x)+1)^3-8/a^3/(tan(1/2*x)+1)^2-20/a^3/(tan(1/2*x)+1)

________________________________________________________________________________________

maxima [B]  time = 0.84, size = 306, normalized size = 3.03 \[ -\frac {\frac {2375 \, \sin \relax (x)}{\cos \relax (x) + 1} + \frac {5347 \, \sin \relax (x)^{2}}{{\left (\cos \relax (x) + 1\right )}^{2}} + \frac {9230 \, \sin \relax (x)^{3}}{{\left (\cos \relax (x) + 1\right )}^{3}} + \frac {12622 \, \sin \relax (x)^{4}}{{\left (\cos \relax (x) + 1\right )}^{4}} + \frac {13340 \, \sin \relax (x)^{5}}{{\left (\cos \relax (x) + 1\right )}^{5}} + \frac {11684 \, \sin \relax (x)^{6}}{{\left (\cos \relax (x) + 1\right )}^{6}} + \frac {8050 \, \sin \relax (x)^{7}}{{\left (\cos \relax (x) + 1\right )}^{7}} + \frac {4370 \, \sin \relax (x)^{8}}{{\left (\cos \relax (x) + 1\right )}^{8}} + \frac {1725 \, \sin \relax (x)^{9}}{{\left (\cos \relax (x) + 1\right )}^{9}} + \frac {345 \, \sin \relax (x)^{10}}{{\left (\cos \relax (x) + 1\right )}^{10}} + 544}{15 \, {\left (a^{3} + \frac {5 \, a^{3} \sin \relax (x)}{\cos \relax (x) + 1} + \frac {13 \, a^{3} \sin \relax (x)^{2}}{{\left (\cos \relax (x) + 1\right )}^{2}} + \frac {25 \, a^{3} \sin \relax (x)^{3}}{{\left (\cos \relax (x) + 1\right )}^{3}} + \frac {38 \, a^{3} \sin \relax (x)^{4}}{{\left (\cos \relax (x) + 1\right )}^{4}} + \frac {46 \, a^{3} \sin \relax (x)^{5}}{{\left (\cos \relax (x) + 1\right )}^{5}} + \frac {46 \, a^{3} \sin \relax (x)^{6}}{{\left (\cos \relax (x) + 1\right )}^{6}} + \frac {38 \, a^{3} \sin \relax (x)^{7}}{{\left (\cos \relax (x) + 1\right )}^{7}} + \frac {25 \, a^{3} \sin \relax (x)^{8}}{{\left (\cos \relax (x) + 1\right )}^{8}} + \frac {13 \, a^{3} \sin \relax (x)^{9}}{{\left (\cos \relax (x) + 1\right )}^{9}} + \frac {5 \, a^{3} \sin \relax (x)^{10}}{{\left (\cos \relax (x) + 1\right )}^{10}} + \frac {a^{3} \sin \relax (x)^{11}}{{\left (\cos \relax (x) + 1\right )}^{11}}\right )}} - \frac {23 \, \arctan \left (\frac {\sin \relax (x)}{\cos \relax (x) + 1}\right )}{a^{3}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(x)^6/(a+a*sin(x))^3,x, algorithm="maxima")

[Out]

-1/15*(2375*sin(x)/(cos(x) + 1) + 5347*sin(x)^2/(cos(x) + 1)^2 + 9230*sin(x)^3/(cos(x) + 1)^3 + 12622*sin(x)^4
/(cos(x) + 1)^4 + 13340*sin(x)^5/(cos(x) + 1)^5 + 11684*sin(x)^6/(cos(x) + 1)^6 + 8050*sin(x)^7/(cos(x) + 1)^7
 + 4370*sin(x)^8/(cos(x) + 1)^8 + 1725*sin(x)^9/(cos(x) + 1)^9 + 345*sin(x)^10/(cos(x) + 1)^10 + 544)/(a^3 + 5
*a^3*sin(x)/(cos(x) + 1) + 13*a^3*sin(x)^2/(cos(x) + 1)^2 + 25*a^3*sin(x)^3/(cos(x) + 1)^3 + 38*a^3*sin(x)^4/(
cos(x) + 1)^4 + 46*a^3*sin(x)^5/(cos(x) + 1)^5 + 46*a^3*sin(x)^6/(cos(x) + 1)^6 + 38*a^3*sin(x)^7/(cos(x) + 1)
^7 + 25*a^3*sin(x)^8/(cos(x) + 1)^8 + 13*a^3*sin(x)^9/(cos(x) + 1)^9 + 5*a^3*sin(x)^10/(cos(x) + 1)^10 + a^3*s
in(x)^11/(cos(x) + 1)^11) - 23*arctan(sin(x)/(cos(x) + 1))/a^3

________________________________________________________________________________________

mupad [B]  time = 7.02, size = 110, normalized size = 1.09 \[ -\frac {23\,x}{2\,a^3}-\frac {23\,{\mathrm {tan}\left (\frac {x}{2}\right )}^{10}+115\,{\mathrm {tan}\left (\frac {x}{2}\right )}^9+\frac {874\,{\mathrm {tan}\left (\frac {x}{2}\right )}^8}{3}+\frac {1610\,{\mathrm {tan}\left (\frac {x}{2}\right )}^7}{3}+\frac {11684\,{\mathrm {tan}\left (\frac {x}{2}\right )}^6}{15}+\frac {2668\,{\mathrm {tan}\left (\frac {x}{2}\right )}^5}{3}+\frac {12622\,{\mathrm {tan}\left (\frac {x}{2}\right )}^4}{15}+\frac {1846\,{\mathrm {tan}\left (\frac {x}{2}\right )}^3}{3}+\frac {5347\,{\mathrm {tan}\left (\frac {x}{2}\right )}^2}{15}+\frac {475\,\mathrm {tan}\left (\frac {x}{2}\right )}{3}+\frac {544}{15}}{a^3\,{\left ({\mathrm {tan}\left (\frac {x}{2}\right )}^2+1\right )}^3\,{\left (\mathrm {tan}\left (\frac {x}{2}\right )+1\right )}^5} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(sin(x)^6/(a + a*sin(x))^3,x)

[Out]

- (23*x)/(2*a^3) - ((475*tan(x/2))/3 + (5347*tan(x/2)^2)/15 + (1846*tan(x/2)^3)/3 + (12622*tan(x/2)^4)/15 + (2
668*tan(x/2)^5)/3 + (11684*tan(x/2)^6)/15 + (1610*tan(x/2)^7)/3 + (874*tan(x/2)^8)/3 + 115*tan(x/2)^9 + 23*tan
(x/2)^10 + 544/15)/(a^3*(tan(x/2)^2 + 1)^3*(tan(x/2) + 1)^5)

________________________________________________________________________________________

sympy [B]  time = 50.54, size = 3288, normalized size = 32.55 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(x)**6/(a+a*sin(x))**3,x)

[Out]

-345*x*tan(x/2)**11/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**
8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*t
an(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 1725*x*tan(x/2)**10/(30*a**3*tan(x/2)**11 +
 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x
/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**
3*tan(x/2) + 30*a**3) - 4485*x*tan(x/2)**9/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**
9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*t
an(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 8625*x*tan(x/2)**8/(
30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2
)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3
*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 13110*x*tan(x/2)**7/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**1
0 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*ta
n(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3)
 - 15870*x*tan(x/2)**6/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2
)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**
3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 15870*x*tan(x/2)**5/(30*a**3*tan(x/2)**1
1 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*ta
n(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*
a**3*tan(x/2) + 30*a**3) - 13110*x*tan(x/2)**4/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/
2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a*
*3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 8625*x*tan(x/2)*
*3/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan
(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*
a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 4485*x*tan(x/2)**2/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)
**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3
*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a*
*3) - 1725*x*tan(x/2)/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)
**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3
*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 345*x/(30*a**3*tan(x/2)**11 + 150*a**3*ta
n(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 138
0*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) +
 30*a**3) - 690*tan(x/2)**10/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*t
an(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 7
50*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 3450*tan(x/2)**9/(30*a**3*tan(x/2)
**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3
*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 1
50*a**3*tan(x/2) + 30*a**3) - 8740*tan(x/2)**8/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/
2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a*
*3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 16100*tan(x/2)**
7/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(
x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a
**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 23368*tan(x/2)**6/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**
10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*t
an(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3
) - 26680*tan(x/2)**5/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)
**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3
*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 25244*tan(x/2)**4/(30*a**3*tan(x/2)**11 +
 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x
/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**
3*tan(x/2) + 30*a**3) - 18460*tan(x/2)**3/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9
 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*ta
n(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 10694*tan(x/2)**2/(30
*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)*
*7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*t
an(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 4750*tan(x/2)/(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*
a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**
5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3) - 1088/
(30*a**3*tan(x/2)**11 + 150*a**3*tan(x/2)**10 + 390*a**3*tan(x/2)**9 + 750*a**3*tan(x/2)**8 + 1140*a**3*tan(x/
2)**7 + 1380*a**3*tan(x/2)**6 + 1380*a**3*tan(x/2)**5 + 1140*a**3*tan(x/2)**4 + 750*a**3*tan(x/2)**3 + 390*a**
3*tan(x/2)**2 + 150*a**3*tan(x/2) + 30*a**3)

________________________________________________________________________________________